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1.  Introduction 


Missile  guidance  integrated  fuzing  (GIF)  technology  seeks  to  elimi¬ 
nate  the  need  for  a  separate  on-board  fuze  by  using  missile  guidance 
parameters  to  calculate  warhead  burst  time.  Guidance  data  can  be 
generated  from  missile  seeker  simulation  code  and  from  missile  flight 
test  data.  Visualization  of  a  missile /target  encounter  is  an  important 
step  in  gaining  an  intuitive  understanding  of  the  interaction  of  guid¬ 
ance  parameters  that  are  used  to  determine  warhead  burst  time. 
Analysis  of  missile  guidance  data  that  resulted  in  suboptimal  fuzing 
performance  is  a  first  step  in  improving  GIF  technology. 


2.  Software  Requirements 

This  report  describes  a  software  program  that  is  a  visualization  tool 
for  guidance  integrated  fuzing  (VTGIF)  data,  as  well  as  the  accompa¬ 
nying  simulation-specific  supporting  software.  The  VTGIF  software 
was  written  in  C++,  Silicon  Graphics  (SGI)  version  4.0,  a  graphical 
user  interface  (GUI)  was  implemented  with  the  Open  Software  Foun¬ 
dation  (OSF)  Motif  toolkit,  version  1.2.4,  and  three-dimensional  view¬ 
ing  and  picking  were  implemented  with  SGFs  Open  Inventor  toolkit, 
version  2.0.1.  VTGIF  requires  an  SGI  graphics  workstation  running 
an  operating  system  of  IRIX  5.3  or  higher.  Guidance  data  were  gener¬ 
ated  for  VTGIF  by  the  Army's  multi-mode  Patriot  simulation 
(MPSIM2)  program.^  An  intermediate  PV~WAVE  program 
(PV~WAVE  is  a  software  package  used  for  analyzing  scientific  data) 
was  required  to  format  and  process  the  guidance  data  generated  by 
MPSIM2  before  it  was  displayed  by  VTGIF.  Figure  1  shows  the  proc¬ 
essing  sequence  for  the  data.  Intermediate  files  were  generated  to  store 
the  data  from  MPSIM2  and  PV~WAVE  code. 


Figure  1. 
Processing  block 
diagram  for 
guidance  data. 


Guidance  Aimpoint 

parameters  calculations 


^MPSIM2  Analyst  Manual,  U.S.  Army  Missile  Command,  Redstone  Arsenal,  Alabama  (August  1993). 
(Confidential) 
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3.  Simulated  Guidance  Data 

The  MPSIM2  program  was  selected  for  generating  guidance  data  be¬ 
cause  it  is  a  mature  GIF  simulation  program  that  has  been  tested  and 
validated.  The  guidance  parameters  generated  by  MPSIM2  that  were 
required  by  VTGIF  are  missile  and  target  position  with  respect  to  the 
launch  location;  missile  yaw,  pitch,  and  roll;  missile  antenna  azimuth 
and  elevation  gimbal  angles  with  respect  to  the  missile  body;  target 
azimuth  and  elevation  orientation  angles  with  respect  to  the  missile 
antenna;  simulated  azimuth  and  elevation  boresight  errors;  simulated 
filtered  (guidance  and  fuzing  algorithms)  azimuth  and  elevation 
boresight  errors;  and  simulated  filtered  range  from  the  missile  to  the 
target.  Also,  missile  flight  time  and  fuzing  range  error  estimates  are 
recorded,  but  are  used  only  by  VTGIF  for  readout  purposes.  A  com¬ 
plete  description  of  these  parameters  can  be  found  in  the  MPSIM2 
Analyst  Manual.^'^ 


4.  Data  Processing 

The  guidance  data  and  monopulse  boresight  errors  generated  by 
MPSIM2  required  processing  before  aimpoint  positions  could  be 
graphically  displayed  on  the  target.  Although  this  processing  could 
have  been  incorporated  into  the  VTGIF  software,  PV~WAVE  was  used 
because  it  simplified  the  required  coding  and  also  provided  a 
hardcopy  graphical  output.  Figure  2  shows  a  simplified,  two- 
dimensional  representation  of  the  typical  coordinate  systems  in  a 
missile/target  encounter.  If  a  third  dimension  were  included  in  the 
diagram,  each  angle  shown  would  require  an  azimuth  and  an  eleva¬ 
tion  angle  component.  MPSIM2  uses  these  general  coordinate  sys¬ 
tems,  but  with  several  unique  coordinate  rotations  required  to  inte¬ 
grate  different  sections  of  the  software.  The  processing  required  to 
account  for  these  unique  rotations  is  standard,  but  it  is  not  discussed 
in  this  report  because  it  is  application  specific.  The  code  is  provided 
in  the  appendix. 

The  aimpoint  of  the  missile  seeker  on  the  target  was  calculated  by 
two  methods.  Simulated  guidance  data  could  be  processed  in  either 
the  earth  {e)  or  the  antenna  {a)  coordinate  system.  Several  assump¬ 
tions  were  made  to  simplify  the  calculations.  I  assumed  that  the  dis¬ 
tance  from  the  center  of  gravity  of  the  missile  to  the  origin  of  the 
antenna  coordinate  system  was  a  constant.  This  assumption  is  rea- 


^MPS1M2  Analyst  Manual,  U.S.  Army  Missile  Command,  Redstone  Arsenal,  Alabama  (August  1993). 
(Confidential) 

^Canh  Ly,  MPSIM  Input  Namelist  Variables,  Army  Research  Laboratory,  unpublished  (January  1995). 
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Figure  2. 
Coordinate 
systems  of 
guidance  data. 


sonable  for  GIF  data,  since  most  or  all  of  the  missile's  fuel  will  have 
been  used  by  the  time  that  fuzing  algorithms  are  engaged.  The 
PV~WAVE  code  also  attempted  to  synchronize  the  simulated  mea¬ 
surements  with  the  known  guidance  parameters.  An  exact  synchro¬ 
nization  would  have  required  an  extensive  analysis  of  the  timing  dia¬ 
grams.  Therefore,  to  simplify  the  analysis,  I  assumed  that  there  was  a 
fixed  processing  delay  that  resulted  in  a  small  timing  error.  I  reduced 
the  resulting  range  error  to  approximately  1  m  by  performing  an 
empirically  calculated  range  correction. 

The  position  and  vector  notations  used  below  both  describe  the  vec¬ 
tors  shown  in  figure  2.  The  only  difference  is  that  the  vector  notation 
describes  the  difference  between  two  guidance  parameters,  and  the 
position  notation  describes  guidance  parameters  that  were  relative  to 
the  initial  launch  position  of  the  Patriot  missile.  In  the  equations  be¬ 
low,  matrices  start  with  the  letter  M  and  are  enclosed  by  brackets, 
and  variables  that  can  describe  both  a  vector  and  an  angular  quantity 
are  denoted  with  a  leading  V  or  6,  respectively.  Both  processing  meth¬ 
ods  require  applying  equations  (1)  to  (3)  to  the  required  missile  and 
target  data  to  calculate  the  line  of  sight  (LOS)  from  the  missile  to  the 
target. 

First,  the  vector  from  the  missile's  center  of  gravity  to  the  missile's 
antenna  must  be  transformed  to  the  earth  coordinate  system: 

D,  =  [M^g]  D,  ,  (1) 

where 


Dj,  =  x’ector  from  the  missile  center  of  gravity  to  the  ori¬ 
gin  of  the  missile  antenna  coordinate  system,  with 
respect  to  the  earth  coordinate  system, 

M^E  =  missile  to  earth  rotation  matrix,  based  upon  missile 
yaw,  pitch,  and  roll,  and 
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=  vector  from  the  missile  center  of  gravity  to  the  ori¬ 
gin  of  the  missile  antenna  coordinate  system,  with 
respect  to  the  missile  coordinate  system. 

To  find  the  position  of  the  missile  in  the  earth  coordinate  system,  we 

can  simply  add  vectors: 

=  ,  (2) 

where 

Mg  =  missile  position  measured  at  the  origin  of  the  an¬ 
tenna  coordinate  system,  with  respect  to  the  earth 
coordinate  system,  and 

Mg^  =  missile  position  measured  at  the  center  of  gravity, 
with  respect  to  the  earth  coordinate  system. 

We  can  now  compute  the  LOS  vector,  Lg,  as  shown  in  figure  2: 

L,  =  T,-M,,  (3) 

where 

Lg  =  LOS  vector  from  the  antenna  to  the  geometric  cen¬ 
ter  of  the  target,  with  respect  to  the  earth  coordinate 
system,  and 

Tg  =  target  position  measured  at  the  geometric  center  of 
the  target,  with  respect  to  the  earth  coordinate 
system. 

For  the  method  based  on  the  antenna  coordinate  system,  equations 

(4)  to  (7)  are  used  to  calculate  the  aimpoint  vector.  First,  the  Foresight 

error  vector  is  calculated  from  monopulse  angular  guidance  data: 

=  I  L  I  [cos  ©g.  cos  ,  sin  ,  sin  dj  ,  (4) 

where 


=  monopulse  Foresight  error  vector  from  the  antenna 
to  the  aimpoint  position,  with  respect  to  the  antenna 
coordinate  system,  labeled  "Foresight"  in  figure  2, 

=  azimuth  angle  monopulse  Foresight  error,  and 

0gy  =  elex  ation  angle  monopulse  Foresight  error. 

The  Foresight  error-  vector  is  then  transformed  to  the  earth  coordi¬ 
nate  system: 

^barCf  ~  ^borca  ' 


(5) 


where 


=  monopulse  boresight  error  vector  from  the  antenna 
to  the  aimpoint  position,  with  respect  to  the  earth 
coordinate  system,  and 

Mae  =  anterma  to  earth  rotation  matrix,  based  upon  azi¬ 
muth  and  elevation  missile  gimbal  angles,  and  mis¬ 
sile  yaw,  pitch,  and  roll. 

Once  transformed  to  the  earth  coordinate  system,  the  boresight  error 
vector  can  be  subtracted  from  the  LOS  vector: 


where 


^ainte  ^borce 


(6) 


^aim  ~  aimpoint  vector  from  the  target  to  the  boresight  vec¬ 
tor,  with  respect  to  the  earth  coordinate  system,  la¬ 
beled  "Aimpoint"  in  figure  2. 

Next,  boresight  error  vector  is  transformed  back  to  the  antenna  coor¬ 
dinate  system: 

=  [Me J  (7) 

where 

^ainia  ~  aimpoint  vector  from  the  target  to  the  boresight  vec¬ 
tor,  with  respect  to  the  antenna  coordinate  system, 
and 

Mea  -  earth  to  antenna  rotation  matrix,  based  upon  azi¬ 
muth  and  elevation  missile  gimbal  angles,  and  mis¬ 
sile  yaw,  pitch,  and  roll. 

Aimpoint  angles  can  now  be  determined  for  the  cross  range  and  el¬ 
evation  components  of 

e„„  =  sin-l(V„-„yiLI),  (8) 

where 

^aim  ~  azimuth  or  elevation  aimpoint  angles,  measured 
with  respect  to  the  earth  coordinate  system. 

For  the  method  based  on  the  earth  coordinate  system,  equations  (9) 
to  (13)  are  required  to  calculate  the  aimpoint  vector.  First,  the  LOS 
vector  is  transformed  to  the  antenna  coordinate  system: 

K  =  [MaeI  h  (9) 
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where 


=  LOS  vector  from  the  antenna  to  the  target,  with  re¬ 
spect  to  the  anterma  coordinate  system,  and 

Ma£  =  antenna  to  earth  rotation  matrix  based  upon  azimuth 
and  elevation  missile  gimbal  angles,  and  missile  yaw, 
pitch,  and  roll. 

The  relative  antenna  beam  angles  are  then  determined  for  the  cross 
range  and  elevation  components  of  LOS: 

e^,  =  sm-i(i,/ILI)  ,  (10) 

where 


=  azimuth  or  elevation  angles  between  the  antenna 
beam  vector  and  the  target  vector,  with  respect  to 
the  antenna  coordinate  system. 

Next,  the  aimpoint  angles  and  vectors  are  calculated  for  the  cross  range 
and  elevation  components  of  : 


n  —  f!  ^  f) 

^aim  ^ant  ^bore  ' 

(11) 

^airria  “  ^  0  / 

(12) 

^ainie  ~  ^ainia 

(13) 

where 

Mra  =  earth  to  antenna  rotation  matrix  based  upon  azimuth 
and  elevation  missile  gimbal  angles,  and  missile  yaw, 
pitch,  and  roll. 


5.  Validation 

The  aimpoint  calculations  were  validated  with  known  guidance  data 
generated  by  MPSIM2.  A  test  target  model  was  constructed  that  con¬ 
sisted  of  a  single  point  scatterer  that  changed  its  position  on  the  tar¬ 
get  as  a  function  of  relative  aspect  angle  with  the  missile  seeker. 
MPSIM2  was  run  in  a  low-noise  environment  with  this  target  model 
and  several  different  flight  paths.  The  results  indicated  that  the  el¬ 
evation  aimpoint  calculations  were  correct,  but  the  azimuth  aimpoint 
calculations  were  dependent  upon  the  relative  sign  of  target  azimuth 
aspect  angle.  The  solution  was  to  change  the  sign  of  the  azimuth  aim- 
point  calculation,  which  worked  satisfactorily.  This  problem  was  stud¬ 
ied  only  briefly  since  the  focus  of  the  effort  was  to  visualize  GIF  data, 
not  to  validate  MPSIM2  (which  is  a  major  task).  Figures  3  and  4  show 
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that  the  aimpoint  position  tracks  the  position  of  a  single  scatterer  on 
the  test  target.  Both  processing  methods  produced  nearly  identical 
results. 


Figure  3.  Validation  of  azimuth  aimpoint  calculations. 


Figure  4.  Validation  of  elevation  aimpoint  calculations. 
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6.  VTGIF  Software 


6.1  Introduction 

GIF  data  were  visualized  with  the  software  program  VTGIF,  which 
allows  the  user  to  examine  the  spatial  relationships  between  target 
aimpoints  and  scattering  centers,  and  between  target  and  missile  po¬ 
sitions  during  a  missile/ target  encounter.  The  software  displays  the 
target,  interceptor,  antenna  beam  pattern,  aimpoint  positions,  and 
target  scattering  centers.  The  display  window  is  typically  centered 
on  the  target,  and  the  target  appears  to  be  stationary  during  most  of 
the  encounter.  When  the  target  orientation  angle  changes  with  re¬ 
spect  to  the  missile  seeker  by  an  amount  greater  than  or  equal  to  the 
angular  resolution  of  the  point  scatter  model,  the  target  will  be  ro¬ 
tated  and  the  corresponding  point  scatterers  will  be  displayed.  Dif¬ 
ferent  target  and  missile  perspectives  can  be  examined  via  a  three- 
dimensional  (3D)  viewer  provided  by  the  Inventor  toolkit.  Individual 
point  scatterer  radar  cross  section  (RCS)  and  position  can  be  displayed 
by  the  user  by  "picking"  them  with  the  mouse.  An  extensive  collec¬ 
tion  of  viewing  formats  and  options  are  available  to  the  user  to  cus¬ 
tomize  the  display.  Figure  5  shows  a  snapshot  of  a  single  frame  from 
the  VTGIF  display. 


Figure  5. 
VTGIF  display 
window. 


ms 
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6.2  Menu  Selections 


The  software  is  controlled  through  a  standard  Motif  menu  interface 
shown  in  the  top  left  comer  in  figure  5.  Different  geometric  target 
models,  point  scatterer  models,  and  simulation  runs  can  be  selected 
through  the  file  menu.  The  geometric  target  model  format  was  a  tri¬ 
angular  mesh  described  by  a  list  of  vertex  coordinates.  The  point  scat¬ 
terer  model  format  was  target  orientation  angle  data,  followed  by  50 
point  scatterers.  This  format  was  repeated  for  each  target  orientation 
angle.  Each  point  scatterer  consisted  of  a  position  in  meters  and  an 
RCS  amplitude  in  meters.  The  simulation  data  consisted  of  a  series  of 
47  guidance  parameters  output  by  the  PV~WAVE  transformation  soft¬ 
ware.  Most  of  the  parameters  were  not  directly  used  by  VTGIR  The 
encoimter  could  be  animated  at  various  speeds  selected  through  the 
animate  menu  interface,  or  stepped  through  by  the  user  frame  by 
frame.  The  target  orientation,  range  to  the  target,  range  error,  and 
total  missile  flight  time  are  displayed  for  each  frame. 

6.3  3D  Viewing 

The  user  can  select  different  3D  viewing  perspectives  using  either  the 
thumb  wheels  shown  in  figure  5  or  the  left  mouse  button.  Different 
viewing  options  can  be  selected  from  a  spring-loaded  menu  activated 
by  the  right  mouse  button  or  by  the  pushbuttons  on  the  right  side  of 
the  window  in  figure  5.  Extensive  on-line  help  is  available  for 
3D  viewing. 


7.  Conclusion 

VTGIF  was  de\'eloped  for  visualizing  a  target /missile  encounter, 
overlaying  target  point  scatterers,  and  displaying  aimpoint  tracking 
positions.  The  software  is  most  useful  for  targets  with  distinct  regions 
where  scatterers  are  concentrated.  For  example,  scatterers  may  be 
concentrated  in  the  nose  section,  along  specific  spots  on  the  body, 
and  on  tail  fins.  Visual  analysis  can  be  used  to  determine  which  sec¬ 
tion  is  being  tracked,  and  what  guidance  parameters  are  responsible 
for  a  target  not  being  adequately  tracked.  The  code  for  VTGIF  is  not 
included  in  this  report  because  it  is  still  under  development.  Future 
software  development  may  include  the  use  of  fragmentation  models 
to  animate  warhead  detonation. 
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Appendix.  Listing  of  PV~WAVE  Code 


;  sanitized  PV~WAVE  program  to  process  MPSIM  GIF  data 
;  antenna  coordinate  system  based-method 
;  I/O  has  been  eliminated  from  this  code 

;  potentially  sensitive  constants  have  been  eliminated  from  this  code 
pro  gif 

cg_m  =  [~cg_len,  0,0]  ;  vector  from  missile  eg  to  antenna 
for  i=l,  n-  delay  do  begin 

r2  =  (data(xm,i-  delay)  -  data(xt,i-  delay)) ^2.  +  (data{hm,i-  delay)  - 

data(ht,i-  delay)  )''2.  +  (data(zm,i-  delay)  -  data(zt,i-  delay))  ^2. 

if  r2  LE  0  then  rg=l 

rg  =  r2^.5  ;  range  from  missile  to  target 

rv  =  [rg,0,0] 

si  =  3.1415926/2.  ;  coordinate  rotation 

si_90  =  [[cos{si),  0,  sin(si)],  $ 

[0,  1,  0],  $ 

[-sin(si) ,  0,  cos(si)]] 

fi  =  pitch_offset 

fi_s  =  [[1,  0,  0],$ 

[0,  cos(fi),  -sin(fi)3,  $ 

[0,  sin(fi),  cos(fi)]] 

si  =  -data {sim, i) ;  yaw 

si_m  =  [[cos (si),  0,  sin{si)],  $ 

[0,  1,  0],  $ 

[-sin(si),  0,  cos(si)]] 

theta  =  -data(thm,i)  ;  pitch 

theta_m  =  [ [cos (theta) ,  -sin(th€ta),  0  ],  $ 

[sin (theta),  cos (theta),  0  ],  $ 

[0,  0,  1]] 
fi  =  -data(fim,i)  ;  roll 
fi_m  =  [ [1,  0,  0]  ,$ 

[0,  cos(fi),  -sin{fi)3,  $ 

[0,  sin(fi),  cos(fi)]) 


theta  =  -data (elgs , i )  ;  antenna  elevation  gimbal  angle 
theta_s  -  [ [cos (theta) ,  -sin(theta) ,  0  3,  $ 

[  sin(theta),  cos (theta), 0  ],  $ 

[0,  0,  133 

si  =  -data (azgs , i )  ;  antenna  azimuth  gimbal  angle 
si_s  =  [[cos(si),  0,  sin(si)3,  $ 

[0,  1,  03,  $ 

[-sin(si),  0,  cos(si)33 


miss_p_e  =  [data(xm,i-  delay) , data (hm, i-  delay) , data (zm, i-  delay) 3  ;  missile  position 
targ_e  =  [data(xt,i-  delay) , data (ht , i-  delay) , data (zt, i-  delay) 3  ;  target  position 

cg_e  =  si_90#si_m#theta_m#f i_m#cg_m 

miss_e  =  cg_e  +  miss_p_e;  ;  corrected  missile  position 

los_e  =  targ_e  -  miss_e  ;  line  of  sight  between  target  and  missile 

;  boresight  vectors 

bore_meas_a= [cos (data (az_meas , i) ) *cos (data (el_meas , i ) ) , sin (data (el_meas, i) ) , - 


sin  (data  (az_ineas,  i) )  ]  *rg 

bore_lak_a  =  [cos (data (az^lak, i) )* cos (data (el^lak, i) ) , sin (data (el^lak, i) 
sin (data (az_lak, i) )] *rg 

bore_fuze_a  =  [cos  (data  {az_fuze,  i) )  ’^cos  (data  (el^fuze ,  i ) )  /  sin  (data  (el_fuze,  i ) ) 
sin  (data  (az_fuze,  i) )]  *rg 

;  convert  from  antenna /seeker  frame  to  earth  frame 

bore_meas_e=si_90#si_m#theta_m#fi_m#fi_s#theta_s#si_s#bore_meas_a 

bore_lak_e=si_90#si_m#theta_m#fi_m#fi_s#theta_s#si_s#bore_lak_a 

bore_fuze_e=si_90#si_m#theta_m#fi_m#fi_s#theta_s#si_s#bore_fuze_a 

aim_dif f_meas_e  =  los_e  -  bore_meas_e 
aim_diff_fuze_e  =  los_e  -  bore_fuze_e 
aim_diff_lak_e  =  los_e  -  bore_lak_e 


sign=l . 

if  data(psif,i)  GT  0  then  sign  =  -1. 

aim_dif f_meas_e (0)  =  aim_dif f_meas_e ( 0) *sign 
aim_diff_fuze_e(0)  =  aim_dif f_fuze_e (0) *sign 
aim„diff_lak_e(0)  =  aim_diff_lak_e (0 ) *sign 

;  convert  aim_diff_e  to  aim_diff_a 

si  =  -3.1415926/2.  ;  coordinate  rotation 

si_90  =  [[cos(si),  0,  sin(si)],  $ 

[0,  1,  0],  $ 

[-sin(si),  0,  cos(si)]] 

si  =:  data  (sim,  i)  ;  yaw 

si_m  =  [[cos(si),  0,  sin(si)],  $ 

[0,  1,  0],  $ 

[-sin(si),  0,  cos(si)]] 

theta  =  data(thm,i)  ;  pitch 

theta„m  =  [ [cos ( theta) ,  -sin(theta),  0  ],  $ 

[sin(theta),  cos (theta),  0  ],  $ 

[0,  0,  1]] 

fi  =  data(fim,i)  ;  roll 

fi^m  =  [  [1,  0,  0) ,  $ 

[0  ,  cos (fi) ,  -sin (fi)  ]  ,  $ 

[0,  sin(fi) ,  cos (fi)  ]  ] 

fi  =  -pitch_of fset 
fi_s=[[l.  0,  0],$ 

[0,  cos(fi),  -sin(fi)],  $ 

(0,  sin(fi) ,  cos(fi) ] ] 

theta  =  data{elgs,i)  ;  antenna  elevation  gimbal  angle 
theta_s  =  ( (cos (theta) ,  -sin(theta),  0  ],  $ 

[  sin (theta),  cos ( theta), 0  ],  $ 

[0,  0,  1]] 

si  =  data(azgs,i)  ;  antenna  azimuth  gimbal  angle 
si_s  =  [(cos{si},  0,  sin(si)],  $ 

[0,  1,  0],  $ 

[ -sin (si )  ,  0,  cos (si)  ]  ] 

;  convert  from  earth  frame  to  antenna/ seeker  frame 

a  im_di  f  f  _meas_a=s  i_s  #  theta_s  #  f  i__s  #  f  i_m#  the  ta„m#s  i_m#  s  i_9  0  #a  im_di  f  f  _meas_e 
aim_diff_lak_a=si_s#theta_s#fi_s#fi_m#theta_m#si_m#si_90#aim_diff_lak_e; 
aim„diff_fuze_a=si_s#theta_s#fi_s#fi_m#theta_m#si_m#si_90#aim_diff_fuze_e; 

;  calculate  aimpoint  angles,  and  assign  values  to  output  array 


16 


data(az_meas, i)  =  asin {aim_dif f_meas_a (2 ) /rg) 
data (el_meas, i)  =  asin (aim_dif f_meas_a (1) /rg) 

data(a2_fuze, i)  =  asin(aim_diff_fuze_a(2) /rg) 
data{el_fuze, i)  =  asin(aim_diff_fuze_a (1) /rg) 

data(az_lak, i)  =  asin(aim_diff_lak_a{2) /rg) 
data(el_lak, i)  =  asin(aim_diff_lak_a(l) /rg) 

return 

end 
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